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Abstract 

We present a general framework for constructing singular solutions of nonlinear evolution equations 
that become singular on a d-dimensional sphere, where d > 1. The asymptotic profile and blowup 
rate of these solutions are the same as those of solutions of the corresponding one-dimensional 
equation that become singular at a point. We provide a detailed numerical investigation of these new 
singular solutions for the following equations: The nonlinear Schrodinger equation iipt(t, x) + Atp + 
\tl)\ 2a i/) = with a > 2, the biharmonic nonlinear Schrodinger equation iipt(t,x.) — A 2 xp + \if;\ 2a 'ip = 
with a > 4, the nonlinear heat equation ipt(t,x) — Aip — \ip\ 2a ip = with a > 0, and the nonlinear 
biharmonic heat equation ipt(t,x.) + A 2 t(j — |-i/>| 2CT, (/j = with a > 0. 



1. Introduction 

In this study, we consider nonlinear evolution equations of the form 

Ut(t,x) = F[u, An, A 2 n,---], x G M. d , d > 1. (1) 

Examples for such equations are the nonlinear Schrodinger equation, the biharmonic nonlinear 
Schrodinger equation, the nonlinear heat equation, and the biharmonic nonlinear heat equation. It 
is well known that these equations admit solutions that become singular at a point. Recently, it 
was discovered that the nonlinear Schrodinger equation with a quintic nonlinearity admits solutions 
that become singular on a d-dimensional sphere [H, 0, IE 0] , see Figure [TJ Following 0] , we refer to 
these solutions as singular standing-ring solutions. 

The main goal of this study is to present a general framework for constructing singular standing- 
ring solutions of nonlinear evolution equations of the form (pQ). In order to understand the basic 
idea, let us assume that equation (pQ) admits a singular standing-ring solution. Then, near the 
singularity, equation (JTJ) reduces to the one-dimensional equation 

Ut(t, r) = F[u, u rr , u rrrr ,--- ], r = |x|. (2) 

Hence, equation ([2]) "should" admit a solution that becomes singular at a point. Conversely, if the 
one-dimensional equation ([2]) admits a solution that becomes singular at a point, then equation (pQ) 
"should" admit a standing-ring singular solution. Moreover, the asymptotic profile and blowup rate 
of the standing-ring solutions of (pQ) "should" be the same as those of the corresponding solution 
of the one-dimensional equation ([2]). 
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Figure 1: Illustration of a two-dimensional singular standing-ring solution at times to < t\ < 
The insets show the radial profile of u. 

The above argument is obviously very informal. Nevertheless, in what follows we will provide 
numerical evidence in support of the relation between standing-ring singular solutions of ([1]) and 
singular solutions of the one-dimensional equation ([2]). 

1.1. Peak-type and ring-type singular solutions of the nonlinear Schrddinger equation (NLS) - re- 
view 

The focusing nonlinear Schrodinger equation (NLS) 

^(t,x) + AV + M 2,T ^ = 0, V(0,x) = Vo(x), (3) 

where x 6 M. d and A = d XlXl + • • • + d XdXd , is one of the canonical nonlinear equations in physics, 
arising in various fields such as nonlinear optics, plasma physics, Bose-Einstein condensates, and 
surface waves. One of the important properties of the NLS is that it admits solutions which become 
singular at a finite time, i.e., 

lim IIV'IIh 1 = °°; < T c < oo. 

t — >Tc 

The NLS is called subcritical if ad < 2. In this case, all solutions exist globally. In contrast, 
solutions of the critical (ad = 2) and supercritical (ad > 2) NLS can become singular in finite time. 
When the initial condition ipo is radially-symmetric, equation ([2D reduces to 

iMt,r) + iPrr + ^—^A + \*P\ 2a *P = 0, VKM = Mr), d>l, (4) 
r 

where r = |x|. Let us denote the location of the maximal amplitude by 

r max (t) = argmax|V>|. 

r 

Singular solutions of (H|) are called 'peak-type 1 when r max (t) = for < t < T c , and 'ring-type' 
when r max (t) > for < t < T c . 
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Until a few years ago, the only known singular NLS solutions were peak-type. In the critical 
case ad = 2, it has been rigorously shown [5( that peak- type solutions are self-similar near the 
singularity, i.e., ip ~ tfin, where 

d-S T fmax ify r,\ „ 

the self-similar profile R{p) attains its global maximum at p = 0, and the blowup rate L{t) is given 
by the loglog law 

In the supercritical case {ad > 2), the rigorous theory is far less developed. However, formal 
calculations and numerical simulations 0] suggest that peak-type solutions of the supercritical 
NLS collapse with the self-similar ips profile, i.e., tfi ~ tps, where 

fe(t,r) = ^SWe'^ 2 , (6a) 
P= F7a ' r max {t) = 0, (6b) 



\S(p)\ attains its global maximum at p = 0, and the blowup rate is a square-root, i.e. 



L(t) ~ k^/T c - t, t -► T c . (6c) 
In the last few years, new singular solutions of the NLS were discovered, which are ring-type 0, 

aaaa. in a 

, we showed that the NLS with d > 1 and | < er < 2 admits singular ring-type 
solutions that collapse with the i/jq profile, i.e., ifi ~ ipQ, where 

^ Q (t, r) = 1 Q (p) e iT+i a ^+i(l-a)±(r-r max (t))2 ^ (?&) 

L ' a (t) 

P = r?T\ ' r max{t) = r L a (t), (7b) 



L\sY ' L(t) 
and 

a = W^Ty (7c) 

The self-similar profile Q attains its global maximum at p = 0. Hence, r max {t) is the ring radius 
and L(t) is the ring width, see Figure [2j 

The ipQ ring solutions can be classified as follows, see Figure [3j 

A. In the subcritical case {ad < 2), all NLS solutions globally exist, hence no singular ring 
solutions exist. 

B. The critical case ad = 2 corresponds to a = 1. Since r max (t) = roL{t), these solutions undergo 
an equal-rate collapse, i.e., the ring radius goes to zero at the same rate as L{t). The blowup 
rate of these critical ring solutions is a square root. 
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Figure 2: Illustration of ring radius r max (i) and width L(t). 



C. The supercritical case 2/d < a < 2 corresponds to < a < 1. Therefore, the ring ra- 
dius r max (i) = roL a (t) decays to zero, but at a slower rate than L(t). The blowup rate of 
these ring solutions is 

L(t) ~ k(T c - tf , (8) 

where p = ^ = 2 +^ } 2) ■ 

D. The supercritical case a = 2 corresponds to a = 0. Since r max (i) = ro, the solution becomes 
singular on the d-dimensional sphere |x| = ro, rather than at a point. The blowup rate of 
these solutions is given by the loglog law . 

E. The case a > 2 was open until now. 

Thus, the tpQ solutions are shrinking rings (i.e., limt~-*T c r max(t) = 0) f° r | < a < 2 (cases B and 
C), and standing rings (i.e., < lim^T c r mas(i) < oo) for a = 2 (case D). 

Singular standing-ring solutions of the NLS 

One of the goals of this paper is to study singular ring- type solutions for a > 2 (case E). The 
most natural guess is that these solutions also blowup with the i/jq profile. Since a < for a > 2, 
see f|7c|> . tpQ should be an expanding ring, i.e., \xmt^>T a r max = oo. In this study we show that 
although such expanding rings do not violate power conservation, tpQ ring solutions cannot exist 
for a > 2. Rather, singular rings solutions of the NLS with a > 2 are standing rings. 

The blowup profile and rate of standing-ring solutions can be obtained using the following 
informal argument. In the ring region of a standing-ring, ip rr ~ and ^ip r ~ — . Therefore, 
the ^^-i^r term in equation becomes negligible compared with ^ rr as t — > T c . Hence, near the 
singularity, equation @ reduces to the one-dimensional NLS0 

i(p(t,x) +(j) xx + \4>\ 2a (p = 0, x = r-r max . 



1 Throughout this paper, we denote the solution of the one-dimensional NLS by <j>, and its spatial variable by x. 
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Figure 3: Classification of singular ring solutions of the NLS as a function of a and d. 

A: subcritical case. B: critical case, with equal-rate collapse 0- C: 2/d < a < 2, shrinking 

rings [2]. D: o" = 2, standing rings 0, El, 0j- E: the case a > 2, which is considered in this 
study. 



Therefore, the blowup profile and blowup-rate of standing-ring solutions of the NLS ((4|) with d > 1 
and cr > 2 are the same as those of singular peak-type solutions of the one-dimensional NLS with 
the same a. Specifically, standing-ring solutions are self-similar in the ring region, i.e., ij) ~ tpp, 
where tpp is, up to a shift in r, the asymptotic peak-type profile of the one-dimensional NLS, i.e., 

ip F (t, r - r max ) = 4> s (t, x), 

and <f>s is given by (I6ap with d = 1. In addition, the blowup rate L{t) of ipp is the same as the 
blowup rate of cps, see (j6c|) . and is equal to a square-root, i.e., 

L(t)~K(a)y/T c -t, t T c . (9) 

Moreover, re(cr) is universal, i.e., it depends on a, but not on the dimension d or the initial condi- 
tion ipQ. 

Numerically, we observe that ipF is an attr actor for a large class of radially-symmetric initial 
conditions, but is unstable with respect to symmetry-breaking perturbations. 

1.3. Singular standing vortex solutions of the NLS 
The two-dimensional NLS 

itpt(t,x,y) + Aip + \ip\ 2a ip = 0, ip(0,x,y) = ip (x,y), A = d xx + d yy , 

admits vortex solutions of the form ifj(t,x,y) = A(t,r)e md , where m G Z. In [3], we presented a 
systematic study of singular vortex solutions. In particular, we showed that there exist singular 



5 



vortex solutions that collapse with the asymptotic profile ipQ ■ e im6 , when ipQ is given by ([7]). The 
blowup rates of these solutions are the same as those of ipQ in the non-vortex case. Therefore, 
the ipQ ■ e l vortex solutions can be classified as follows: 

A. In the subcritical case (a < 1), all NLS solutions globally exist, hence no singular vortex 
solutions exist. 

B. In the critical case a = 1, these solutions undergo an equal-rate collapse at a square root 
blowup rate. 

C. The supercritical case 1 < a < 2 corresponds to < a < 1. In this case, the ring radius decays 
to zero at a slower rate than L(t) and the blowup rate is given by ([8]) where p = = §. 

D. The supercritical case a = 2 corresponds to a = 0. Therefore, the solution becomes singular 
on a circle. The blowup rate is given by the loglog law ([5]). 

E. The case a > 2 was open until now. 

In this study, we show numerically that there exist singular standing-vortex solutions of the 
two-dimensional NLS with a > 2 (case E) with the asymptotic profile e m8 tpF- Moreover, the 
blowup rate of these singular standing vortices is the same as that of the standing-ring solutions in 
the non- vortex case, i.e., is given by Q. Therefore, these results extend the ones obtained in the 
non- vortex case. 

1.4- Singular solutions of the biharmonic nonlinear Schrddinger equation 

Let us consider the focusing biharmonic nonlinear Schrodinger equation (BNLS) equation 

# i (t,r)-A 2 V + |V'| 2CT V' = 0, (10) 

where A 2 is the radial biharmonic operator. Here, singularly formation is defined as lim | IV 1 1 1 #2 = 

00. In the subcritical case ad < 4, all BNLS solutions exists globally 0]. Numerical simulations [1,0] 
indicate that in the critical case a = 4/d and the supercritical case a > 4/ci, the BNLS admits 
singular solutions. At present, however, there is no rigorous proof that the BNLS admits singular 
solutions, whether peak-type or ring-type. 

Peak-type singular solutions of the BNLS (|10p were recently studied numerically in [9]. The 
blowup rate of these solutions is slightly faster than p = 1/4 in the critical case (1/4+loglog?), and 
is equal to p = 1/4 in the supercritical case. 

The BNLS (110p also admits ring-type singular solutions for 4/d < a < 4 [0j]. These solutions 
are of the form ip ~ ipQ B , where 

\^Qb\ = TTT^T^Qb (p) , (Ha) 



LV2*(t) 



L(Jt) 

and 



,(t)=r Q L aB (t), (lib) 



aB = ~H — TV ( llc ) 
a(d — 1) 

The ipQ B solutions can be classified as follows (see Figure 0J: 

A. In the subcritical case (ad < 4), all BNLS solutions globally exist, hence no collapsing ring 
solutions exist. 
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B. The critical case ad = A corresponds to = 1. These solutions undergo an equal-rate 
collapse. The blowup rate of these critical ring solutions is given by ([8]) with p = 1/4. 

C. The supercritical case A/d < a < A corresponds to < as < 1. Therefore, the ring ra- 
dius r max (i) = roLg(t) decays to zero, but at a slower rate than L(t). The blowup rate of 
these ring solutions is given by (jHJ) with p = 1/(3 + olb) = &(d — l)/(4 + 3a d — Aa). 

D. The case a = A corresponds to «b = 0. S ince ?" max (t) = tq, the solution is a singular standing 
ring. The blowup rate is close top = 1/4 and is conjectured to be 1/4 with a loglog correction. 

E. The case a > A was open until now. 
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Figure 4: Classification of singular ring solutions of the BNLS as a function of a and d. A: 
subcritical case. B: critical case, with equal-rate collapse [9j. C: A/d < a < 4, shrinking 
rings . D: er = 4, standing rings with the critical ID profile j!) . E: the case a > A, which is 
considered in this study. 



Thus, up to the change a — ► 2a, this classification is completely analogous to that of singular ring 
solutions of the NLS (see Figure [3]). In this work we show numerically that this analogy carries 
through to the regime a > A. Thus, the BNLS with a > A and d > 1 admits singular standing-ring 
solutions. Near the standing-ring peak, equation (|10p reduces to the one-dimensional BNLS 

i(f>t(t,x) - <\> xxxx + |0| 2<T = 0. 

Therefore, the blowup profile and blowup-rate of standing ring solutions of the BNLS (|10p with 
d > 1 and a > 4 are the same as those of collapsing peak solutions of the one-dimensional BNLS 
with the same value of a. 

Thus, the results for ring solutions of the BNLS with a > A are completely analogous, up to 
the change a — * 2a, to those for singular standing-ring solutions of the NLS with a > 2. 
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1.5. Singular solutions of the nonlinear heat equation and the biharmonic nonlinear heat equation 
The (i-dimensional nonlinear heat equation (NLHE) 



ut(t, r) — Au 



\u\ 2a u 



0, 



a > 0, d > 1, 



(12) 



and the d-dimensional nonlinear biharmonic heat equation (NLBHE) 

u t (t,r) + A 2 u - \u\ 2a u = 0, a > 0, d > 1. 



(13) 



admit singular solutions for any cr > 1CJ, lUJ] . To the best of our knowledge, until now all known 
singular solutions of (|12|) and of (|13|) collapse at a point (or at a finite number of points [12]). In 
this study, we provide numerical evidence that the NLHE (|12h and the NLBHE fjl3[> admit singular 
standing-ring solutions. The blowup profile and blowup rate of these solutions are the same as 
those of singular peak-type solutions of the corresponding one-dimensional equation. 



1.6. Critical exponents of singular ring solutions 

In Figure [5A we plot the blowup rate parameter p of singular ring solutions of the NLS, see ([8]). 
As a increases from 2/d to 2—, p increases monotonically from 1/2 to 1—. When a = 2, the blowup 
rate is given by the loglog law ([5]), i.e., p = 1/2 with a loglog correction. Finally, p = 1/2 for a > 2. 
Since 

lim p = 1, lim p = — , 

the blowup rate has a discontinuity at cr = 2. Surprisingly, the blowup rate is not monotonically- 
increasing with a. For example, a ring solution of the NLS with a = 1.8 blows up faster than a 
ring solution of the NLS with a = 2.2. 

The above results show that the critical exponent of singular ring solutions of the NLS is a = 2: 
The blowup rate is discontinuous at a = 2, and the blowup dynamics changes from a shrinking- 
ring (a < 2) to a standing-ring (cr > 2), see Figure [5p. We can understand why a = 2 is a critical 
exponent using the following argument. Standing-ring solutions are 'equivalent' to singular peak 
solutions of the one-dimensional NLS with the same nonlinearity exponent a. Since a = 2 is the 
critical exponent for singularity formation in the one-dimensional NLS, it is also a critical exponent 
for standing-ring blowup. 

An analogous picture exists for the BNLS. For A/d < a < 4, the blowup-rate p of the BNLS 
ring solutions increases monotonically in a from 1/4 to (1/3) — , at a = 4, p = 1/4 possibly with a 
loglog correction, and p = 1/4 for a > 4, see Figure [5jB. Therefore, the blowup rate is discontinuous 
at cr = 4. In addition, the blowup dynamics change at cr = 4 from a shrinking-ring (cr < 4) to 
a standing-ring (cr > 4), see Figure [5p. Hence, the critical exponent of standing-ring solutions of 
the BNLS is a = 4, precisely because it is the critical exponent for singularity formation in the 
one-dimensional BNLS. 

In the case of the NLHE and BNLHE equations, there is no critical exponent of singular 
ring solutions. Indeed, these equations admit standing-ring solution for any cr > 0, precisely 
because there is no critical exponent for singularity formation in the corresponding one-dimensional 
equations. 
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Figure 5: A: Blowup rate of singular ring solutions of the NLS. The blowup rate increases 
monotonically from p = 1/2 at a = 2/d to p = 1— at a = 2— . For a = 2 (full circle) p = 1/2 
(with a loglog correction) and for a > 2, p = 1/2. B: Blowup rate of singular ring solutions of 
the BNLS. The blowup rate increases monotonically from p = 1/4 at a = 4/d to p = (1/3)— 
at a = 4—. For a = 4 (full circle) p = 1/4 (with a loglog correction?) and for a > 4, p = 1/4. 
C: The shrinkage parameter a, defined by the relation r max ~ roL a of singular ring solutions 
of the NLS. For 2/d < a < 2, a decreases monotonically from 1 to 0+ (shrinking rings). 
For a > 2, a = (standing rings). D: The shrinkage parameter as of singular ring solutions 
of the BNLS. For A/d < a < 4, as decreases monotonically from 1 to 0+ (shrinking rings). 
For a > 4, ob = (standing rings). 
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1.7. Paper outline 

The paper is organized as follows. In Section [2] we review the theory of singular peak-type 
solutions of the supercritical NLS, and conduct a numerical study of the one-dimensional case. In 
Section \3. II we prove that standing-ring blowup can only occur for a > 2, and show that the blowup 
profile and blowup-rate of singular standing ring solutions of the supercritical NLS with a > 2 
and d > 1 are the same as those of peak-type solutions of the one-dimensional NLS equation. 
In Section 13.21 we confirm these results numerically. We then show numerically that the singular 
standing-ring profile ipF is an attractor for radially-symmetric initial conditions (Section 13.30 . but 
it is unstable with respect to symmetry-breaking perturbations (Section 13. 4p . In Section [J] we show 
analytically and numerically that expanding ipQ ring solutions do not exist for a > 2. Section [5] 
extends the results to singular vortex solutions. In Section [6] we study singular peak-type solutions 
of the one-dimensional supercritical BNLS. In Section [7] we show that singular ring solutions of 
the supercritical BNLS with a > 4 are standing-rings, whose blowup profile and blowup-rate are 
the same as those of peak-type solutions of the one-dimensional BNLS. In Section [8] we show that 
singular standing-ring solutions of the nonlinear heat equation exist for any a > 0, and that their 
blowup profile and blowup-rate are the same as those of peak-type solutions of the one-dimensional 
NLHE. In Section [9] we show that singular standing-ring solutions of the nonlinear biharmonic heat 
equation exist for any a > 0, and that their blowup profile and blowup-rate are the same as those 
of peak-type solutions of the one-dimensional BNLHE. The numerical methods used in this study 
are briefly described in Section PTOl 



2. Singular peak-type solutions of the one-dimensional supercritical NLS 

2.1. Theory review 

Let us consider the one-dimensional supercritical NLS 

i<Mt,x) + <fe x + H 2<T = O, er>2. (14) 

In contrast to the extensive theory on singularity formation in the critical NLS, much less is 
known about the supercritical case. Previous numerical simulations and formal calculations (see, 
e.g., [6, Chapter 7] and the references therein) suggested that peak-type singular solutions of the 
supercritical NLS (fT4l) collapse with a self-similar asymptotic profile <jf>s, .i.e., (j) ~ cps, where 



The blowup rate L(t) of these solutions is a square root, i.e., 



L(t) ~ Ky /T c - t, t -► T c , (16) 
where k > 0. In addition, the self-similar profile S is the solution of 

In general, solutions of (|17p are complex- valued, and depend on the parameter k and on the initial 
condition S(0). Solutions of (|17p whose amplitude |»S| is monotonically-decreasing in £, and which 
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Figure 6: The parameters k$ and 5q of the admissible solutions of equation (|17|) . as function 
of a. 



have a zero Hamiltonian, are called admissible solutions [g]. For each a, equation (fT7|) has a unique 
admissible solution (up to a multiplication by a constant phase e ia ). This solution is attained for 
specific real values of k and 5(0), which we denote as 

k = K S (a), 5(0) = 50(a). (18) 

Moreover, numerical simulations and formal calculations suggest that: 

1. The self-similar profile of singular peak- type solutions of the NLS (1141) is an admissible solution 
of (HZD- 

2. The constant k of the blowup rate (|16p is universal (i.e., is independent of the initial condi- 
tion V'o)) an d is equal to Ks'(cr). 

2.2. Simulations 

To the best of our knowledge, the theory of supercritical peak-type collapse which is presented 
in Section \2. 11 was tested numerically only for d > 2. Since this theory is not rigorous, and since we 
will make use of these results in Sections [3] and HJ we now confirm numerically the above theoretical 
predictions for the one-dimensional supercritical NLS (|14p . 

We verified numerically that for each <r, there exists a unique admissible solution of (|17p . The 
corresponding values of ks and So as a function of a are shown in Figure [6l For example, for a = 3, 
the parameters of the admissible solution of equation (|17j) are 



k s (ct = 3) ** 1.664, 5 (ct = 3) » 1.155. (19) 



We now solve the one-dimensional NLS (|14p with a = 3 and the Gaussian initial condi- 
tion cfio (t = 0, x) = 2e~ 2x . We first show that the NLS solution collapses with the self-similar 
profile 05, see (]15p . To do that, we rescale the solution according to 



Waled (t, x) = L*(t) 



x 



L(t) 
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Figure 7: Solution of the one-dimensional NLS (|14p with a = 3 and the initial condi- 
tion ifjQ = 2e~ 2x . A: Rescaled solution, see (f20j) . at the focusing levels 1/L = 10 4 (solid) 
and 1/L = 10 8 (dashed); dotted curve is the admissible profile S(£;cr = 3), all three curves 
are indistinguishable. B: L as a function of (T c — t) on a logarithmic scale. Dotted curve is the 
fitted curve 1.713 • (T c — i) - 5007 . The two curves are indistinguishable. C: LLt as a function 
of 1/L. 



The rescaled profiles at focusing levels of 1/L = 10 4 and 1/L = 10 8 are indistinguishable, see Fig- 
ure [7]A, indicating that the solution is indeed self-similar while focusing over 4 orders of magnitude. 
Moreover, the rescaled profiles are in perfect fit with the admissible S(£; a = 3) profile. 
Next, we consider the blowup rate of (p. To do that, we first assume that 

L(t) ~ k(T c - if, (21) 

and find the best fitting k and p, see Figure 03. In this case k ~ 1.713 and p ~ 0.5007, indicating 
that the blowup rate is square-root or slightly faster. In order to check whether L is slightly faster 
than a square root, we compute the limit limLLj, see [7(j. Recall that for the square-root blowup 

rate (fi"6|) . 

K 2 

lim LL t = — 7T < 0, 

t^T c 2 

while for a faster-than-a-square root blowup rate LLt goes to zero. Since lim LLt = —1.384, see 
Figure [7)0, the blowup rate of 4> is square-root (with no loglog correction), i.e., 

L(t) ~ K \ l ° WUp - latC y / T~^t, K ^wup-rate _ y^OM « 1.664. (22) 

In particular, there is an excellent match (to 4 digits) between K ^ wu P~ ratc ~ i.g64 extracted from 
the blowup rate of (p, see ([22]) and the parameter K5(cj = 3) of the admissible S((,;a = 3) profile, 
see (USD. 



3. Singular standing-ring solutions of the supercritical NLS 

Let us consider singular solutions of the supercritical NLS 

iipt{t,r) + ip Tr + ^—^ip r + \i)\ 2(T ip = {), d > 1, ad>2. (23) 
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In this Section we show that equation (I23p has singular standing-ring solutions for a > 2. Since 
the case a = 2 was already studied in 0, El, 0] , we mainly focus on the case a > 2. 

3.1. Analysis 

The following Lemma shows that standing-ring blowup can only occur for a > 2: 

Lemma 1. Let ip be a standing-ring singular solution of the NLS (f2"3"j) . i.e., ip ~ ipp f or \ r ~ ?"maxl < 
p c - Lit), where 



\ipF(t,r)\ 



and 



< lim r max (i) < oo. 
t^T c 



lim L(t) = 0, 

t^T c 



(24a) 
(24b) 



Then, a > 2. 

Proof. The power of the collapsing part t/>_f is 



Vila 



F 



r - r n 



L 



r d-l dr 



>r=rmax—pc-L(t) 

L- 2/a r \F(p)\ 2 (Lp + r max ) d - 1 (Ld /9 ) 

J p=-p c 

rpc 



rX-2/a(j.\ . d-1 
^i/j 'max 



P=-pc 



Since ||^ F || 2 < M 1 1 2 
Let 



('olio < °°> then L 1 2 /°" has to be bounded as L — ► 0, hence <r > 2. □ 



-Pcoiiapse = liminf lim / 



|^|V _1 dr 

be the amount of power that collapses into the standing-ring singularity. We say that a singular 
standing-ring solution ip undergoes a strong collapse if -P co iiapse > 0, and a weak collapse if -P co iiapse = 
0. 

Corollary 2. Under the conditions of LemmaUl ipp undergoes a strong collapse when a = 2, and 
a weak collapse when a > 2. 



Proof. This follows directly from the proof of Lemma [TJ 



□ 



Let us further consider singular standing-ring solutions of the NLS (|23p . Following the analysis 
in [1, section 3.2], in the ring region of a standing-ring solution, i.e., for r — r max = O(L), 



[lp r 



L 2 (ty 



d-l 



'"max 



Therefore, the ^zr-ip r term in equation (I23j) becomes negligible compared with ip rr as f — > T c . 
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Hence, near the singularity, equation ([23]) reduces to the one-dimensional supercritical NLS {HJ, 

i.e., 

ip(t, r) ~ <f>(t,x = r - r max (t)) , 

where <f> is a peak-type solution of the one-dimensional NLS (1141) . 

Therefore, we predicted in [2] that the blowup dynamics of standing ring solutions of the 
NLS (123p with d > 1 and cr = 2 is the same as the blowup dynamics of collapsing peak solutions 
of the one-dimensional critical NLS with a = 2, as was indeed confirmed numerically in [2] and 
analytically in [H, 0]. Similarly, we now predict that the blowup dynamics of standing-ring solutions 
of the NLS (|23p with d > 1 and a > 2 is the same as the blowup dynamics of collapsing peak 
solutions of the supercritical one-dimensional NLS (|14p with the same nonlinearity exponent a: 

Conjecture 3. Let if) be a singular standing-ring solution of the NLS equation (f23j) with d > 1 
and a > 2. Then, 

1. The solution is self-similar in the ring region, i.e., \ip\ ~ \iPf\ for r — r max = 0(L), where tpp 
is given by (I24al) . 

2. The self-similar profile ipp is given by 

ip F (t, r) = <f>s(t,x = r - r max (t)) , (25) 

where 4>s(t,x), see \15\). is the asymptotic peak-type profile of the one- dimensional NLS 
with the same a. In particular, F = S is the admissible solution of equation (fT7|) with 

k = k s (g), So = S (cr). 

3. The blowup rate of ip is a square root, i.e., 

L(t) ~ K S (a)VT c -t, t — > T c , (26) 
where Ks(o~),is the parameter k = ks of the self-similar profile S, see (|18p . 

Conjecture implies that the parameter n of the blowup-rate (|26|) of tp is equal to the parameter 
k of the blowup rate (|16p of (p. In particular, k depends on the nonlinearity exponent a, but is 
independent of the dimension d and of the initial condition ?/V 

3.2. Simulations 

We solve the NLS (|23p with d = 2 and a = 3 for the initial condition 

V>o = 2e- 2 ^- 5 ) 2 , (27) 

and observe that the solution blows-up with a ring profile. In Figure [8] we plot the ring radius 

?"max(i) = argmaxlV'l 
r 

as a function of the focusing factor 1/L(i), as the solution blows up over 10 orders of magnitude. 
Since lim^r c r max (t) = 5.0011, the ring is standing and is not shrinking or expanding. 
We now test Conjecture [3] numerically item by item. 
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Figure 8: Ring radius r max {t) as a function of the focusing level 1/L(t) for the solution of the 
NLS (f23|) with d = 2 and a = 3, and the initial condition ([27]) . 



1. In Figure 0A. we plot the rescaled solution 

T Tmax (^) 



V'rescaled — L"(t) 



L(t) 



m = iM^) , (28) 



at 1/L = 10 4 and 1/L = 10 8 , and observe that the two lines are indistinguishable. Therefore, 
we conclude that the standing-ring solutions blowup with the self-similar ipp profile (|24aj) . 

2. To verify that the self-similar blowup profile ipp is> U P to a shift in r, the asymptotic blowup 
peak-profile 4>s of the one-dimensional NLS (|14p . we superimpose in Figure [9]A. the self-similar 
profile of the solution of the one-dimensional NLS (|14|) from Figure [7]A and the admissible 
solution S(x,a = 3), and observe that, indeed, the four curves are indistinguishable. 

3. Figure [9(3 shows that 

L(t) ~ 1.714- (T c -t) a5009 . 
Therefore, the blowup rate is a square-root or slightly faster. Figure[9p shows that lim LLt ~ 

T c >t 

— 1.385, indicating that the blowup rate is square-root (with no loglog correction), i.e., 

L(t) ~ ^ wup - rate v / TW, K^ wup " rate = V2 • 1.385 « 1.664. (29) 

Thus, there is an excellent match between the parameter k = ks(<t = 3) ~ 1.664 of the 
admissible S profile, see (fl~9|) . the value of K2^ wup rate ~ 1.664 extracted from the blowup rate 
of the solution of the two-dimensional NLS, and the value of K ^ wu P~ ratc ~ 1.664 extracted 
from the blowup rate of the solution of the one-dimensional NLS, see (|22p . 



3.3. Robustness ofipp and universality of k 

The initial condition (|27p in Figures [8] and [9] is different from the asymptotic profile ipp- Since 
the solution ip blows up with the asymptotic profile ipp, this indicates that tpp is an attractor. 
The initial condition (|27p . however, is already ring-shaped. Therefore, we now show that initial 
conditions which are not ring-shaped can also blowup with the V'f profile. 
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Figure 9: NLS solution of Figure [8l A: Rescaled solution according to (|28p at focusing levels 
l/L = 10 4 (solid) and l/L = 10 8 (dashed), dotted curve is the asymptotic profile 5(£, a = 3), 
and the dashed curve is the rescaled solution of the one-dimensional NLS at l/L = 10 8 , taken 
from Figure OA.. All four curves are indistinguishable. B: L as a function of (T c — t) on a 
logarithmic scale. Dotted curve is the fitted curve 1.709(T C — i) - 5007 . C: LLt as a function of 
l/L. 



In 13, ll4|, we developed a nonlinear Geometrical Optics (NGO) method which showed that 
high-power super-Gaussian initial conditions evolve into a ring profile. To see this, in Figure [TOl we 
solve the NLS (J3j) with d = 2 and a = 3, and the super-Gaussian initial condition tpo(r) = 2e~ r , 
and observe that the NLS solution, indeed, evolves into a ring. Since \\m.r max (t) = 0.33, see 

Figure [TTTA.. this singular solution is a standing ring. Therefore, we see that initial conditions which 
are not rings can also blowup with the ipp standing ring profile. 

We now consider the blowup rate of the above solution, since lim LL t = —1.384 , see FigurefTTB. 

this implies that 

L(t) ~ 4 l ° WUP ~ VatC V / T~^t, K blowup-rate _ ^J84 = 1.664. 

This value of At2^ wup_rate identifies with the one obtained for the ring- type initial condition (|27p . 
see (j29|) . We thus see that the parameter k of the blowup rate ([26]) is, indeed, independent of the 
initial condition. 

Remark 3.1. A different type of initial condition that blows-up with the ipF profile and with the 

i p blowup— rate • • n , ■ 1—5 — 75 

same value of k 2D is given in bection \4-ty 

3.4- Instability with respect to symmetry-breaking perturbations 

In Section [3.31 we saw that the standing-ring asymptotic profile ipF is an attractor for a large 
class of radially-symmetric initial conditions. In general, NLS solutions with a ring structure are 
stable under radial perturbation, but unstable under symmetry-breaking perturbations 0, S Sj- 
We now show that i^f is also unstable with respect to symmetry-breaking perturbations. To see 
that, let us consider the two-dimensional NLS 



iMt, r, 0) + Ar + + -$>w + l^l 2 "^ = 0, (30) 
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Figure 10: Solution of the NLS (|23p with d = 2 and a = 3 and the initial condition -00 = 2 e 
at A: t = 0. B: t = 0.0105. C: t = 0.0188. D: t = 0.0198. 
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Figure 11: NLS solution of Figure [TU1 A: Location of maximum r max (t) as a function of the 
focusing level 1/L(t) B: LL^ as a function of the focusing level 1/L(t). 
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with the initial condition 

Mr,e) = f(r)(l + eh(9)). 

We chose f(r) so that when e = 0, the solution blows up with the t/jp profile at r = r max . We now 
consider the case < e < 1. Since for a standing ring the ^ip r term becomes negligible compared 
with ip rr , see Section [37TT equation ([30]) can be approximated in the ring-peak region (r « r max ) 
with 

ii/H + Ar + -J— Vtoo + IV^I 2ct -V> = 0. 

This is the two-dimensional focusing NLS. Therefore, the solution will localize at local maximum 
points in the (r, 6) plane, thereby breaking the radial symmetry. 

To see this numerically, we solve the two-dimensional NLS (I30p with a = 3 and with the initial 
condition 

^(r,0)=2e- 2 ( r - 5 ) 2 [l + e 2 e-(f) 2 ], e=± = [-tt,tt]. (31) 

This initial condition is the standing-ring initial condition (|27p . with an O(0. 01) small bump 
at 6 = 0, see Figure PT21A. As predicted, as the solution self-focuses, it localizes around the small 
initial bump at = (see Figure [T2B and!12C). resulting in breakup of radial symmetry. 



4. Existence of non-standing ring solutions for a > 2 ? 

Lemma [T] does not exclude the possibility that there exist non-standing rings for a > 2. The 
main reason that this question arises is as follows. In [2] we discovered ring solutions of the 
supercritical NLS for d > 1 and | < a < 2 of the form ip ~ ^q, see ([7]). Therefore, it is natural 
to attempt to extrapolate these results to the regime a > 2. Since a < when a > 2, the ring 
radius r max (t) goes to infinity as t — > T c , hence ipQ is an expanding-ring profile for a > 2, if such 
a solution exists. Note that although the ring is expanding to an infinite radius, the power of the 
collapsing part ipQ remains bounded, as 

\^ Q \ 2 r d - l dr ~ rt 1 / \Q\ 2 dp, t -> T c . 

=r max -p c -L(t) Jp=-p c 

Therefore, these expanding rings, if they exist, do not violate power conservation. 

In [3] we solved the NLS (|23l) with a = 2.1 > 2 and <i = 2 and the ring initial condi- 
tion = v / 3\/sech(2(r — 5)) . The solution turned out to be a singular standing ring, rather 
than an expanding one. Moreover, the blowup profile was different from ipQ. In retrospect, this 
NLS solution was a standing-ring with the ipF profile, see Section [3l Nevertheless, this still leaves 
open the question of whether there exist expanding-ring ^jq solutions for a > 2. 

4-1. Analysis 

We now prove that singular ring solutions with the ipQ profile do not exist for a > 2. 
Lemma 4. When a > 2, there are no singular NLS solution such that tp ~ ipQ, see (JY]), and 

L(t)~ K(T c -t) p , L t ~ P K(T c -t) p -\ L tt ~ p (p-l)K(T c -t) p - 2 . (32) 
Proof. The result shall follow directly from Lemmas and EJ □ 
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Figure 12: Solution of the NLS (I30p with a = 3 and with the initial condition (J3T]). A: t = 0. 
B: t = 0.01400, C: t = 0.01472. Top: Surface plot. Bottom: Amplitude along the ring 
peak \ip(t,r max (t),9)\ as function of 6. 
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Lemma 5. Under the assumptions of Lemma^ p < 1. 



Proof. We first recall that, as shown by Merle [151 ] . for every singular solution ifi of the supercritical 
NLS 

T 

(T c -t)||W||!<#<°°- (33) 

o 

To find the limiting behavior of ||V^(i)||| a s t — > T c , note that by the conservation of the Hamil- 
tonian 

1 

ff + l 

In addition, 



l|v^l!~— MiKS. t-r c . (34) 



/ \Q{p)\ 2<T+2 {Lp + r L«) d - l Ldp~ ^ + l aL l {t) J \Q\ 2a ' 



d-l 

2a+2 ||,/, \\2a+2 

2a+2 ~ 2a 



l dp, (35) 



where in the last equality we used the value of a given by ([Tc"]) . and, in particular, that a < 0. 
Therefore, by ([32]). (|3l|). P|l 

»v*li ~ MISS. ~ ~ pr^, * - r„ 

Hence, the bound ([33]) implies that p < 1. □ 
Lemma 6. Under the conditions of Lemma^ p > 1. 

Proof. Substitution of ^q, see ([7]), into the NLS ([23]) gives the following ODE for Q, 

Qpp(p) + L+^a Qp ~Q + \Q\ 2a Q ~ [ A ^)P 2 + »r B(t)p + ar 2 C(t)} Q + iD{t)Q = 0, (36a) 



where 



A{t) = h 3 L tt , B(t) = \L 2+a L tt - 2(1 - a)L 1+a L 2 , 
C(t) = -L l+2a L tt - (1 - a)L 2a L 2 , D(t) = ^1 ^—p. 



Since Q depends only on p, each of the time-dependent terms of ([36]) should go to a constant 
as t — > T c . In particular, C{t) should go the constant as t — > T c . Under assumption (j32[) . 



C(t) ~ c c (T c - t ) 2a P+ 2 P- 2 , c = K 2+2a p 



(37) 



Since lim^^ L(t) = 0, then p > 0, see ([32]) . In addition, for a > 2, a < 0, see ([7c]) . This implies 
that c < and in particular c / 0. Since C(i) should goes to a constant as t — ► T c then, by ([37]) . 
2ap + 2p — 2 > . Therefore, 

1 — v 

< a < 0. 

p 

Hence, p > 1. □ 
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Figure 13: Solution of the NLS (I23p with d = 2 and a = 3 and the initial condition (|38f) . A: 
Ring radius r max {t) as a function of the focusing level 1/L(t). B: LL^ as a function of 1/L. 



4-2. Simulations 

The result of Lemma H] that expanding-ring singular solutions with the profile ipQ do not exist, 
is based on formal arguments rather than on a rigorous analysis. Therefore, we now provide a 
numerical support for this result. To do that, we solve the NLS (I23p with d = 2 and a = 3 and 
with the expanding ring profile initial condition 

iP = ip Q (t = 0) = (1 + a)^sech^(cj(r - i )) e -^ 2 - i ( 1 - Q )('- 10 ) 2 , (38a) 

where 

9 — 3 1 

a= , \ =--. (38b) 
3(2-1) 3 y ' 

If a tpQ solution indeed exists, then ip would be a singular ring solution whose radius goes to 
infinity. In Figure [T3lA we plot the ring radius r max (t) as a function of the focusing factor 1/L, 
as the solution blows up over 10 orders of magnitude. Initially, the ring radius, indeed, expands 
from r max (0) = 10 to 

Tmaxif) ~ 12.11. This expansion is due to the defocusing (expanding) phase 
term e~ tar of the initial condition. However, the ring stops to expand when 1/L ~ 20, and becomes 
a singular standing ring with radius r max (T c ) ~ 12.11. Since the initial condition was an expanding 
ring, this simulation provides a strong support to the result of Lemma [H 

We now consider the blowup rate of the above solution, FigurefT3B shows that lim LLt = —1.384, 

t >T C 

implying that 

L(t) ~ K ^° wup ~ iatc v / T~^t, />£ wup - rate w V2 • 1.384 = 1.664. 

This value of K2i) WUp-rate identifies with the one obtained for a ipp collapse, see (|29p . Therefore, 
this simulations provides an additional support to the robustness of ipp and to the universality of k 
(see Section l3T3j) . 

5. Singular standing vortex solutions of the NLS (er > 2) 

We now consider vortex solutions of the two-dimensional NLS 

iMt,x,y) + A^ + \ij\ 2a ^ = 0, 1>(P,x t y)=ifo(x,y), A = d xx + d yy , (39) 
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i.e., solutions of the form 

ip(t, r, 9) = A(t, r)e me , m G Z, (40) 



where r = \J x 1 + y 2 and 9 = tan 1 {x/y). 

In H we proved that if the initial condition is a radially-symmetric vortex, then the solution 
remains a vortex: 

Lemma 7. Let tp be a solution of the NLS Ii39\) with the initial condition Tpo(r,6) = Ao(r)e im . 
Then, ip(t,r,9) = A(t,r)e l , where A(t, r) is the solution of 

iA t (t,r)+A rr + ^A r -^ T A + \A\ 2a A = 0, A(0,r) = A {r). (41) 

Note that the phase singularity at r = implies that A(r = 0) = 0. Hence, all vortex solutions 
are ring-type solutions. Specifically, all the singular solutions of (|4"Tj) are ring-typed and not peak- 
typed. 

In [3J we showed by formal calculations and numerical simulations that equation (|39p admits 
singular shrinking- vortex solutions for 1 < a < 2, and singular standing- vortex solutions for a = 2. 
Moreover, we showed that the blowup rate and profile of the standing-vortex solutions is the same 
as in the two-dimensional non- vortex case. We now show that this is also true for a > 2, namely, 
that the analysis conducted in Section 13.11 for non-vortex standing-ring collapse, applies also for 
singular standing-vortex solutions. 

5.1. Analysis 

Lemma 8. Let ijj be a singular standing-ring vortex solution of the NLS (j39[) . i.e., ijj ~ i^Fit, r)e ime , 
where i/jp is given by (j24a|) . Then, a > 2. 

Proof. The proof is identical to the proof of Lemma [TJ Indeed, the proof of Lemma Q] relies only 
on \i[>\, hence is not affected by the phase term e im8 . □ 

We now show that the blowup dynamics of standing vortex solutions is the same as the blowup 
dynamics of collapsing solutions of the one-dimensional NLS (|14p . Indeed, in the ring region of a 
standing vortex solution, 



[Ar 



[A] 



(d-l)[A] 

?max(7c) - L(t) 



1 

m . 
-^A 



m 2 [A] 



Therefore, as t — > T c , both the ^r-A r and ^-A terms in equation (123j) become negligible compared 
with A rr . 

Hence, as in the non-vortex case, near the singularity, equation (|23p reduces to the one- 
dimensional NLS (|14|) . i.e., 

A{t,r)~<p(t,x = r-r max (t)), (42) 



where <j) is a peak-type solution of the one-dimensional NLS (|14p . Therefore, we expect that the 
blowup dynamics of standing-vortex solutions of the NLS (|23p with d = 2 and a > 2 to be the same 
as the blowup dynamics of collapsing peak solution of the one-dimensional NLS equation (|14p with 
the same nonlinearity exponent a: 
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Figure 14: Ring radius r max (t) as a function of the focusing level 1/L(t) for the solution of 
the two-dimensional NLS (|23p with m = 1, a = 3 and the initial condition (|43p . 



Conjecture 9. Lei r, 0) = j4(t, r)e* &e a singular standing-vortex solution of the NLS equa- 
tion (|39p wii/i a > 2. Then, tp blows up with the asymptotic self- similar profile 

1> ~ e im0 -i; F (t,r), 

where ipp is given by (|25p . In addition, items 1-3 of Conjecture^ hold. 

5.2. Simulations 

We solve equation (|4ip with m = 1 and <r = 3, with the initial condition 

^0 = 2 tanh(4r 2 )e- 2 ( r " 5 ) 2 . (43) 

In Figure [T4l we plot the ring radius r max {t) as a function of the focusing factor 1/L(t), as the 
solution blows up over 10 orders of magnitude. Since limt_,.r c r max (t) = 5.0011, the vortex is 
standing. 

We now test Conjecture [9] numerically item by item. 

1. In Figure [15] we plot the rescaled solution, see equation (i28j) . at 1/L = 10 4 and 1/L = 10 8 , 
and observe that, indeed, the standing ring solution undergoes a self-similar collapse with the 
profile (p5|). 

2. To verify that the self-similar collapse profile is, up to a shift in r and multiplication by e imS , 
the asymptotic collapse profile efts of the one-dimensional NLS (|14p . we superimpose the 
rescaled solution of the one-dimensional NLS (|14p from Figure [7K, as well as the admissible 
solution S(x,a = 3), onto the rescaled solutions of Figure [T51A and observe that, indeed, the 
four curves are indistinguishable. 

3. Figure [T5B shows that 

L(t) ~ 1.701 • (T c -t) a50068 . 
Therefore, the blowup rate is square root or slightly faster. Figure [T51G shows that 

lim LL t w -1.384, 
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Figure 15: Solution of equation (|4ip with a = 3, m = 1, and the initial condition (|43|) . 
A: Rescaled solution according to (l28l) at focusing levels 1/L = 10 4 (solid) and 1/L = 10 8 
(dashed), dotted curve is the asymptotic profile 5, and the dashed curve is the rescaled so- 
lution of the one-dimensional NLS at 1/L = 10 8 , taken from Figure [7]A. All four curves are 
indistinguishable. B: L as a function of (T c — t) on a logarithmic scale. Dotted curve is the 
fitted curve 1.71(T C — t) ' 5007 . The two curves are indistinguishable. C: LLt as a function of 
1/L. 



indicating that the blowup rate is square-root, i.e., 

L(f) ~ 4»-Z; t T ^-ZZT = ^L384 « 1.664. 

In addition, there is an excellent match between the parameter k = ks(ct = 3) ~ 1.664 of the 
admissible S profile, see (fT9|) . and the value of ^D-vortcx* ~ 1-664 extracted from the blowup 
rate of the two-dimensional vortex solution. 



6. Singular peak-type solutions of the one-dimensional supercritical BNLS 

In Section[2]we reviewed the theory of singular peak-type solutions of the one-dimensional NLS. 
In this section, we present the analogous findings for the one-dimensional BNLS. We will make use 
of these results in the study of singular standing-ring solutions of the BNLS in Section [3 

6.1. Analysis 

Let us consider the one-dimensional supercritical focusing BNLS 

i(f> t (t, x) - <p xxxx + \(f>\ 2r7 <f) = 0, a > 4. (44) 

At present, there is no theory for singular peak-type solutions of equation (|44[) . A recent numerical 
study [s| suggests that peak-type singular solutions of the supercritical BNLS (f4"4"l) collapse with a 
self-similar asymptotic profile (f>B, -i.e., <j) ~ 4>Bi where 

" uk) me ^ s = Wr T(t) = L^k- (45) 
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Figure 16: Solution of the one-dimensional BNLS (I44p with a = 6 and the Gaussian initial 
condition ([4"8|h A: Rescaled solution according to ((15]) at focusing levels 1/L = 10 4 (solid) 
and 1/L = 10 8 (dashed). The two curves are indistinguishable. B: L as a function of (T c — t) 
on a logarithmic scale. Dotted curve is the fitted curve 1.020(T C — i) - 25017 . C: L 3 L^ as a 
function of 1/L. 



The blowup rate L(i) of these solutions is a quartic root, i.e., 

L(t)~n B \/T c -t, t^T c , (46) 
where kb > 0. In addition, the self-similar profile L?(£) is the solution of 

("1 + ^4) B(0 + ^4eB? - + \B\ 2u B = 0. (47) 

Note that the NLS analogue of the B(£) profile is not the S(£) profile of equation (fT7j) . but rather 

.it 2 

• e 4 4r x . The analogue of the quadratic phase term in BNLS theory is at present unknown. 
In general, symmetric solutions of f|4Tj) are complex-valued, and depend on the parameter kb 
and on the initial conditions B(0) and B"(0). We conjecture that, in analogy with the NLS, the 
following holds: 

Conjecture 10. 1. The nonlinear fourth-order ODE (|47p has a unique 'admissible solution' 
with k b = k b {u), B{0) = B {a), and B"{0) = B%{a). 

2. This admissible solution is the self-similar profile of the asymptotic peak-type blowup profile 
4>B, see (1431). 

3. The value of kb of the blowup rate P6]) is equal to kb(o~) of the admissible B profile. 
6.2. Simulations 

We solve the one-dimensional BNLS (I44p with a = 6 and the Gaussian initial condition 

<j)(t = 0,x) = 1.6 -e~ x2 . (48) 

We first show that the BNLS solution blows up with the self-similar profile 4>b, see (I45p . To do 
that, we rescale the solution according to 

descaled (t,x) = L 2 ^(^ ( -^-] , L(t) = \\cf>\\^ 2 . (49) 
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The rescaled profiles at focusing levels of l/L = 10 4 and 1/L = 10 8 are indistinguishable, see 
Figure [TBIA. indicating that the solution is indeed self-similar while focusing over 4 orders of mag- 
nitude. 

Next, we consider the blowup rate of the collapsing solution, see Figure [TBB. To do that, we 
first assume that 

L(t)~K B (T c -t) p , (50) 

and find the best fitting kb and p. In this case kb ~ 1.020 and p as 0.25017, indicating that the 
blowup-rate is close to a quartic root. To verify that the blowup rate is indeed p = 1/4, we compute 



the limit limL Lj. Note that for the quartic-root blowup rate (1461) 

t^T c 

lim L 3 L t = — ^ < 0, 
t^T c 4 

while for a faster-than-a-quartic root blowup rate L s Lt — ► 0. Since lim Tc ^ t L 3 L t a -0.2898, sec 
Figure [TBIC. the blowup rate is a quartic-root (with no loglog correction), i.e., 

L (t) ~ «S°m P_rate ^TW, 4% up " rate « ^4 • 0.2898 a 1.0376. 

7. Singular standing-ring solutions of the supercritical BNLS 

In Section [3] we analyzed singular standing-ring solutions of the NLS with a > 2. In this section, 
we derive the analogous results for the biharmonic NLS with a > 4. 

7.1. Analysis 

Let us consider singular solutions of the focusing supercritical BNLS 

r) - + |-0| 2ct V = 0, ad>4, d>l, (51) 



(d-l)(d-3)_ , (d-l)(d-3)_ 2 , 2(d-l)_3 



where 

A 2 = — — — ^ " ; d r + - — ^ " y d r 2 + ^ — ^d A r +d$ (52) 

is the radial biharmonic operator. The following Lemma, which is the BNLS analogue of Lemma [H 
shows that standing-ring collapse can only occur for a > 4: 



Lemma 11. Let ip be a self-similar standing-ring singular solution of the BNLS (|10p . i.e., ip ~ Vs; 
where 

\iJ>B(t,r)\ = L 2/<r( t \ \ B( <P) I' P= T L(t) — ' ^ 53a ' ) 



and 



< lim r max (i) < oo. (53b) 
t >T C 



Then, a > 4. 

Proof. Integration gives H^sll^ = C (^ 1 4 ^ cr ), and so the proof of Lemma [T] holds for a > 4. □ 
Corollary 12. undergoes a strong collapse when a = 4, and a weak collapse when a > 4. 
Proof. This follows directly from the proof of Lemma [TT1 □ 
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Let us further consider standing-ring solutions of the BNLS (|5ip . In this case, in the ring region 
of a standing ring solution, the terms of the biharmonic operator, see (I52p . behave as 



^flfy =0[L~ k ), fc = 0,...,4. 

Therefore, ~ <9 4 V>- Hence, near the singularity, equation (f5Tj) reduces to the one-dimensional 
BNLS ([Ml), i.e., 

^(t,r)~0(t,a; = r-r max (t)), (54) 

where 4> is a peak- type solution of the one-dimensional BNLS (|44p . Therefore, we predicted in 0] 
and also confirmed numerically that the blowup dynamics of standing ring solutions of the NLS (|5ip 
with d > 1 and a = 4 is the same as the blowup dynamics of singular peak solutions of the 
one-dimensional BNLS (|5ip with a = 4. Similarly, we now predict that the blowup dynamics of 
standing ring solutions of the BNLS (I5ip with <i > 1 and cr > 4 is the same as the blowup dynamics of 
collapsing peak solutions of the one-dimensional BNLS (|44p with the same nonlinearity exponent a: 

Conjecture 13. Let ip be a singular standing-ring solution of the BNLS (|51|) with d > 1 and a > 4. 
Then, 

1. T/ie solution is self-similar in the ring region, i.e., ip ~ tps for r — r max = 0(L), where |i/>b| 



is gruen fry (j53a[) . 

2. T/ie self-similar profile ips is given by 

i>B{t, r) = (j) B (t,x = r - r max (t)) , (55) 

where (f>B(t,x), see ( L/5| ) ; is the asymptotic profile of the one- dimensional BNLS \44\ ) with the 
same a. 

3. The blowup rate is a quartic root, i.e., 



L{t) ~ n B (a) ?jT c -t, t^ T c , (56) 
where Kb(o~) > is the value of Kb of the admissible B profile. 
Conjecture 1131 implies, in particular, that the parameter k-q of the blowup-rate of tjj, see ()56[) . 



is the same as the parameter k-q of the blowup rate of (j), see (|46p . This value depends on the 
nonlinearity exponent a, but is independent of the dimension d and of the initial condition tpQ. 

7. 2. Simulations 

We solve the BNLS (|5ip with d = 2 and a = 6 with the initial condition Vo( r ) = 1-6 • e — ( r ~ 5 ) 2 . 
In Figure 1171 we plot the ring radius function of the focusing factor 1/L(t), as the 

solution blows up over 8 orders of magnitude. Since lim^j^ r max (i) = 4.992, the ring is standing 
and is not shrinking or expanding. Note that the initial condition is different from the asymptotic 
profile ipB, suggesting that standing-ring collapse is (radially) stable. 

We next test each item of Conjecture 1131 numerically: 

1. In Figure PTBlA. we plot the rescaled solution 

VWaicd = L2/a ^( r ~^ {t) )> Ht) = UU /2 , r max (t) = argmaxH, (57) 

at 1/L = 10 4 and 1/L = 10 8 . The two curves are indistinguishable, showing that standing 
rings undergo a self-similar collapse with the self-similar ipg profile (|45p . 
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Figure 17: Ring radius r max (t) as a function of the focusing level 1/L(t) for the solution of 
the BNLS (ID]) with d = 2, a = 6 and the initial condition ip (r) = 1.6 • e~( r ~ 5 ) 2 . 




2. To verify that the self-similar collapse profile ipB is, up to a shift in r, equal to the asymptotic 
collapse profile (f>B of the one-dimensional BNLS (|44j) . we superimpose the rescaled solution of 
the one-dimensional BNLS (|44j) from Figure [TolA. onto the rescaled solutions of Figure [T81A. 
and observe that, indeed, the curves are indistinguishable. 

3. Figure [T8B shows that 

L(t) ~ 1.020(T C - t) - 25017 . 
Therefore, the blowup rate is quartic root or slightly faster. Figure [THD shows that lim L 3 Lt ~ 
—0.2894, indicating that the blowup rate is quartic-root (with no loglog correction), i.e., 

L (t) ~ K^ up - ratc ^IH, ^lowup-rate w ^4 • 0.2894 « 1.0373. 



As predicted, there is an excellent match between the value of k^>°2]j P rat ° ~ 1-0373 extracted 

from the two-dimensional BNLS solution, and the value of k^°^ p ratc ~ 1.0376 extracted 
from the one-dimensional BNLS solution, see Section 16.21 



8. Singular standing-ring solutions of the nonlinear heat equation 

The d-dimensional radially-symmetric nonlinear heat equation (NLHE) 

u t (t,r) - Au-\u\ 2a u = 0, a>0, d > 1, (58) 

where u is real and A = d rr + ^r-d r , admits singular solutions for any a > [hJ]. To the best of 
our knowledge, until now all known singular solutions of (|58p collapsed at a point (or at a finite 
number of points [12]). We now show that the NLHE admits also singular standing-ring solutions. 
The blowup profile and blowup rate of these solutions are the same as those of singular peak-type 
solutions of the one-dimensional NLHE with the same a. 
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Figure 18: Solution of the two-dimensional BNLS (jlOl) with a = 6 and the initial con- 
dition ^o(x) = 1.6 • e~( r_5 ) . A: Rescaled solution according to at focusing levels 
1/L = 10 4 (solid) and 1/L = 10 8 (dashed). The dashed curve is the rescaled solution of 
the one-dimensional BNLS at 1/L = 10 8 , taken from Figure [TBI All three curves are indis- 
tinguishable. B: L as a function of (T c — t) on a logarithmic scale. Dotted curve is the fitted 
curve 1.020(T C — i) - 2502 . The two curves are indistinguishable. C: L 3 Lt as a function of 1/L. 



8.1. Peak-type solutions of the one- dimensional NLHE (review) 
one-dimensional NLHeH 

v t (t, x) - v xx - \v\ 2a v = 0, a > 0, 

admits singular solutions that collapse with the self-similar peak-type profile 

1 1 



where 



v s (t,x) 



A(t) = ^/2a{T c - t) 



\«{t)(l+?)-° 



L(ty 



L(t) 



2(2 + i)c/;.-/)|]uc/;.-/)|. 



(59) 



(60a) 



(60b) 



see |16f | . Note that unlike the one-dimensional NLS, the one-dimensional NLHE admits singular 
solutions for any a > 0. 

8.2. Analysis 

Let us consider singular standing-ring solutions of the NLHE (|58p . Since Au ~ u rr in the ring 
region, near the singularity equation (|58|) reduces to the one-dimensional NLHE (|59|) . Therefore, 
we conjecture that singular standing-ring solutions of the NLHE (|58p exist for any a > 0, and that 
the blowup profile and blowup rate of these solutions are the same as those of singular peak-type 
solutions of the one-dimensional NLHE with the same a. 



2 Throughout this paper, we denote the solution of the one-dimensional NLHE by v, and its spatial variable by x. 
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Figure 19: Solution of the NLHE $M§ with d = 2, a = 3 and the initial condition (|62l) . 
A: i = 0. B: t = 0.002683. C: t = 0.002686. 




Figure 20: Solution of Figure [T9l A: Ring radius r max (i) as a function of 1/L(t). B: Solution u 
at 1/L = 10 8 (solid curve). Dashed curve is the asymptotic us profile (|64p . Both curves are 
rescaled according to (|63p . The two curves are indistinguishable. 



Conjecture 14. Zei u(t,r) be a singular standing-ring solution of the NLHE (|58|) . Then, the 
solution is self-similar in the ring region, i.e., u ^ us for r — r max = O(L), where 

u s (t,r)=v s (t,x = r-r max (t)), (61) 

and vs is given by equation (|60p. 

5.5. Simulations 

We solve the NLHE (|58p with d = 2 and cr = 3 and the initial condition 

uo(r) = 2e" 2(r - 5)2 . (62) 

The solution blows up with a ring profile, see Figure [T9l Since limt^T c r max {t) = 4.9994 > 0, see 
Figure [20lA the ring is standing. 
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In order to show that the solution blows up with the self-similar us profile (|6ip . we rescale the 
solution according to 

"reseated (*> r ) = A-0)« (~J^y^j ' ( 63 ) 

where X(t) and L(t) are given by (|60bp and T c is extracted from the numerical simulation. The 
rescaled profile at 1/L = 10 8 is in perfect agreement with the rescaled us profile 

1 

"S.rescaled = f, (64) 

(1 + p 2 )^ 

see Figure I2UB. Therefore, the numerical results provide a strong support to Conjecture [HJ 



9. Singular standing-ring solutions of the nonlinear biharmonic heat equation 

The (i-dimensional radially-symmetric biharmonic nonlinear heat (BNLHE) equation 

u t (t,r) + A 2 u- \u\ 2(7 u = 0, a>0, d > 1, (65) 
where u is real and A 2 is the radial biharmonic operator (|52p . admits singular solutions for any a > 







To the best of our knowledge, all known singular solutions of the BNLHE (|65p collapse 
at a point. We now show that the BNLHE admits singular standing-ring solutions. The blowup 
profile and blowup rate of these solutions are the same as those of singular peak-type solutions of 
the one-dimensional BNLHE with the same a. 



9.1. Peak-type solutions of the one- dimensional BNLHE (review) 
The one-dimensional BNLHE equation 

v t (t, x) + v xxxx - \v\ 2a v = 0, a > 0, (66) 

admits singular peak-type solutions which collapse with the self-similar peak-type profile 

v ^ x ) = j2j^ B (0, L(t) = K BH ^/T~~t, C = ^y ; (67) 

see fill ]. The self-similar profile -£>(£) is not known explicitly. Unlike the one-dimensional BNLS, 
the one-dimensional BNLHE admits singular solutions for any a > 0. 

9.2. Analysis 

Let us consider singular standing-ring solutions of the BNLHE (|65p . Since A 2 u ~ u rrrr in the 
ring region, near the singularity the BNLHE (|65p reduces to the one-dimensional BNLHE (|66p . 
We therefore conjecture that standing-ring solutions of equation (|65p exist, and that their blowup 
profile and blowup rate is the same as those of singular peak solutions of equation ([66]) : 

Conjecture 15. Let u(t,r) be a singular standing-ring solution of the BNLHE (|65p. Then, the 
solution is self-similar in the ring region, i.e.,u ~ ub for r — r max = 0{L), where 

UB(t,r) = v B (t,x = r - r max (t)) , (68) 

and vb is given by equation (I67p . 
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1/L 

Figure 21: Standing ring solution of the BNLHE (165p with initial condition (|62h . A: ring 
radius r max (i) as a function of 1/L(t) B: solution u rescaled according to ([571) . at focusing 
factors of 1/L = 10 4 (solid curve) and 1/L = 10 8 (dashed curve). The rescaled peak-type 
solution of the one-dimensional NLH (|66|) is given by the dotted curve. All three curves are 
indistinguishable. 
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9.3. Simulations 

We solve the BNLHE (|65p with d = 2 and <r = 3 and the initial condition (|62p . The solution 
blows up with a standing-ring profile, see Figure I2TTA In FigureETjB we plot the solution, rescaled ac- 
cording to (|57p . at focusing levels of 1/L = 10 4 and 1/L = 10 8 . The two curves are indistinguishable, 
indicating that the solution is indeed self-similar. Next, we want to show that the self-similar blowup 
profile is given by -£>(£), the self-similar profile of peak- type solutions of the one-dimensional NLHE, 
see (I67p . To do that, we compute the solution of the one-dimensional BNLHE (|66p with a = 3 
and u(t = 0, x) = 2e~ x , and superimpose its profile at 1/L = 10 8 in Figure l21~B. The three rescaled 
solutions are indistinguishable, indicating that standing-ring solutions of the BNLHE blowup with 
the self-similar profile of peak-type solutions of the one-dimensional BNLHE. In addition, we have 
from the numerical simulations that limt—>T c L 3 Lt ~ —1.2108 when d = 2, where L(t) := [|u||oc > > 
and limt—>Tc L^Lt ~ —1.2108 when d = 1, where L(t) := Hulloo^ 2 - Therefore, the blowup rate of 
the standing-ring solution of the two-dimensional BNLHE is equal, up to 5 significant digits, to the 
blowup rate of the singular peak-type solution of the one-dimensional BNLHE, and is given by 

L{t) ~ k bh {/T c - t, k bh ~ ^4 • 1.2108 w 1.4835. 

Thus, the numerical results provide a strong support for Conjecture [151 



10. Numerical Methods 



10.1. Admissible solutions of (|17p 

The admissible solution S of (|17p was calculated using the shooting method of Budd, Chen and 
Russel [17j, Section 3.1]. In this method, one searches in the two-parameter space (S(0),f c ) for the 
parameters such that the solution of (|17p satisfies the admissible solution condition 



lim F(£;/ c ,S(0)) = 0, F(£;/ c ,S(0)) 



s(0 
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10.2. Solution of the NLS, BNLS, NLHE and BNLHE 

In this study, we computed singular solutions of the NLS ©, the BNLS (JED), the NLHE (|58|) 
and the BNLHE (|65p . Theses solutions become highly- localized, so that the spatial scale-difference 
between the singular region r — r max = O(L) and the exterior regions can be as large as 0(1/L) ~ 
10 10 . In order to resolve the solution at both the singular and non-singular regions, we use an 
adaptive grid. 

We generate the adaptive grids u sing the Static Grid Redistribution (SGR) method, which was 
first introduced by Ren and Wang [18|], and was later simplified and improved by Gavish and 
Ditkowsky 191] . Using this approach, the solution is allowed to propagate (self-focus) until it starts 
to become under-resolved. At this stage, a new grid, with the same number of grid-points, is 
generated using De'Boors 'equidistribution principle', see 3, 3] for details. 



The method in [19(] also allows control of the portion of grid points that migrate into the singular 
region, preventing under-resolution at the exterior regions. In [9J], we further extended the approach 
to prevent under-resolution in the transition region 0{L) <C r — r max <C O(l). 

On the sequence of grids, the equations are solved using a Predictor-Corrector Crank-Nicholson 
scheme. 
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